Multiplicity Fluctuations in Au+Au Collisions at RHIC 
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In the PHOBOS experiment, charged particles are measured in almost the full solid angle. This 
enables the study of fluctuations and correlations in the particle production over a very wide kine- 
matic range. In this paper, we show results of a direct search for fluctuations identified by an 
unusual shape of the pseudorapidity distribution. In addition, we use analysis of correlations of 
the multiplicity in similar pseudorapidity bins, placed symmetrically in the forward and backward 
hemispheres, to test the hypothesis of production of particles in clusters. 



1 Introduction 

In the study of heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC), a new type of strongly 
interacting matter was discovered [1]. During a phase transition, rapid changes of the system's properties are 
possible and unusual fluctuations may be observed[2]. Early predictions of signals of a quark-giuon plasma (QGP) 
suggested the possibility of a large increase of total multiplicity or correlated particle production in limited phase 
space. Although already the first results from RHIC have shown only a moderate increase of charged particle 
multiplicity with collision energy, it is still possible that less apparent but significant fluctuations can be found. 
A detailed analysis of fluctuations is also a source of information on the particle production mechanism in the 
early stage of the collision and in the later hadronization, both of which are important for understanding the 
properties of this new state of matter. 

In our analysis, Au+Au collisions at ^/s NN =200 GeV, the highest energy for heavy nuclei available at RHIC, 
are used. In the first study, the most central (3%) collisions are examined in order to identify events with an 
unusual angular distribution of produced particles. The second study analyses multiplicity fluctuations in forward 
and backward pseudorapidity bins and shows an enhancement with respect to purely statistical expectations. 
This may be explained by production of particles in clusters, which is modeled and compared to the experimental 
data. 

The PHOBOS detector[3] consists of several elements located mostly near the beam pipe. Most important for 
the present analysis is the multiplicity detector, which has a uniquely large geometrical acceptance for measuring 
produced particles. The hits from charged particles emitted at \r)\ < 5.4 are registered in a single layer of 
silicon sensors, arranged as an octagonal tube in the center of the apparatus (for < 3) and mounted in 
three pairs of rings perpendicular to the beam direction for detection of particles emitted at smaller angles. 
We measure the ionization energy deposited by the particle in the silicon and thus can estimate the number of 




Figure 1: Rejection of pileup events with large signal in the ZDC counters (above the line) for all central events 
(left) and events with unusual dN/drj shape (right). 



particles hitting the same pad by dividing the signal by the value expected for a single particle. The multiplicity 
can be reconstructed using this estimation ("analog" method) or by counting the hits and applying occupancy 
corrections ("hit counting" method). The two methods give consistent results and were successfully used to 
measure angular distribution of charged particles in the full acceptance of the PHOBOS detector[4, 5]. 



2 Search for unusual events 

In the previous study [6] of the most central events (3%), an excess of events with large multiplicities was found, 
whose rate was however correlated with luminosity. The fraction of events with a dN/drj distribution deviating 
from the mean shape also increased with the event rate. The most probable explanation of such results is the 
registration of two Au+Au collisions as one event (pileup). In this paper, we apply a more efficient rejection of 
pileup events and focus the analysis on the dN/drj distribution. For each of over 1,900,000 most central events we 
determine the raw number of particles as a function of pseudorapidity, M{rj), using the "hit counting" method 
without occupancy corrections. We also calculate the mean value {M(jj)) and variance crj\ {rj) in narrow bins of 
vertex position (to minimize acceptance fluctuations). The deviation of the 77 distribution for an event is then 
measured by: 

2 1 v (M(r,) - S{M( V )) f 

where A/^bias denotes the number of bins in 77 in which M(rj) is defined and S is the fitted scaling factor accounting 
for varying total multiplicity. We define as unusual the events with X^df a va l ue which should not be exceeded 
in our sample. In the previous study, pileup events were identified by comparison of signals in the silicon sensors 
(long signal integration time) with signals in the scintillator trigger counters (short signal integration time) . This 
allows to eliminate merging of Au+Au collisions from different bunch crossings. In this paper, we also compare 
signals in the silicon sensors with the signal left by spectator neutrons in the Zero Degree Calorimeters (ZDC). 
For the most central Au+Au collisions, the ZDC signal should be small, thus the events above the lines in Fig. 1 
are rejected, as most probably the central collision is accompanied by a peripheral one. This cut removes over 
90% of the events with x^ DF >3 and, together with the previous cut, reduces the tail of large X^- DF significantly 
(Fig. 2). Properties of the remaining unusual events will be investigated in a future study, but we know already 
that the fraction of central events with large deviations of their dN/drj shape is smaller than 10~ 5 . 



3 Forward-backward fluctuations 

Another aspect of multiplicity fluctuations is studied by comparing the multiplicity measured in the forward and 
backward hemispheres, reconstructed event-by-event in pseudorapidity bins of varying size and position. The 
multiplicity N F reconstructed in a bin in the forward rj region is compared to the multiplicity N B in a similar 
bin placed symmetrically at negative rj using a correlation parameter C: 
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Figure 2: Distribution of X^df ^ or a ^ cen t ra l events (histogram) and these accepted by the pileup cuts (shaded 
area) . 



which is calculated as a function of position of the center of the positive rj bin or of its width, A77. In the 
comparison of bins covering similar kinematical regions the expected mean value of C is zero. If, in addition, 
multiplicity fluctuates only statistically according to the binomial partitioning of N F +N B , we obtain the variance 
cr 2 (C) = 1. Values of cr 2 (C) different from 1 would indicate the presence of nonstatistical fluctuations, possibly 
due to correlations between primary particles. 

In this analysis, a simplified version of the "analog" method based on the energy deposited by charged 
particles in the octagonal multiplicity detector is used to determine 7V F and N B . Using detector simulations, the 
value of <J 2 {C) is corrected for modifications due to secondary production, Landau energy loss fluctuations, and 
acceptance effects [7]. The reconstructed values of c 2 (C) are larger than 1, as is expected in the case of particle 
production in two steps. For a simple example, we can assume that initially produced objects (clusters) decay 
into k particles, all of which are registered in the same 77 bin. If the primary clusters are produced randomly 
(giving cr 2 (C c i us tcr)=l), the variance of C for final particles grows to er 2 (C) = k. If the distribution of cluster 
size, k, has a mean value of (k) and variance (7k, the relevant value is k e g = (k) + (T 2 /(fc). In a more realistic 
case, particles originating from the same cluster are distributed in pseudorapidity, and thus usually are not all 
contained in the same rj bin used in the calculation of C. In addition, if the forward and backward bins are close, 
particles from a single cluster may be emitted into opposite bins. These two effects decrease the value of c 2 (C) 
in a bin-width and bin-position dependent way. 



4 Simple Cluster Models 

The hypothesis of cluster production in Au+Au collisions can be tested using simple models generating clusters 
which then decay into particles. Any such model should reproduce the general properties of Au+Au events: 
the experimental pseudorapidity distribution[5] and the transverse momentum distribution [8], However, perfect 
agreement is not necessary, as we have found that modifications of the dN/dt] shape or the (p t ) value cause only 
small changes of the impact of clusters. 

The Isotropic Cluster Model (ICM) assumes that clusters have negligible transverse momentum and decay 
into particles emitted isotropically in the clusters' rest frame. In this frame, the pseudorapidity distribution of 
particles has a shape similar to a Gaussian function, but with longer tails. The cr(rj) of isotropic clusters is about 
1.0, and is moderately modified by the Lorentz transformation of the final particles to the laboratory frame, 
especially for the largest \r]\ values (it is about 15-20% smaller at |ry|=3). In this model, the multiplicity of a 
cluster can be freely selected. In order to obtain a continuous spectrum of k c s values, we generated k from a 
Poisson distribution with an appropriate mean value. 

Exactly isotropic decay of each cluster is an extreme case, in which (7(77) is the largest. Such angular symmetry 
is not preserved in decays of resonances, which are the most natural candidates for clusters. The momentum 
of an emitted particle, for example a pion in a decay R{M{) — ► i?(Ma) + 7r, depends mostly on the mass 
difference AM = Mi — M2. Only the values AM < 650 MeV lead to momenta of pions smaller or consistent 
with the experimental mean p t . In our Resonances Cascade Model (RCM) we select a chain of decays starting 
from relatively light particles and thus small AM in decays: 

A(1900) — ► A(1600) + 7r — > A(1232) + 2tt — > N + 3ir (3) 



o rrci Q 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 0.5 1 1.5 2 2.5 3 



Ar| m 

Figure 3: Comparison of reconstructed cr 2 (C) with the predictions from two cluster models: the dependence on 
the width of the r\ bin (left) and on the position of the bin center (right). 

In the later analysis we assume that one of the final particles is neutral, thus the size of the cluster is always 3. 
In order to reproduce even approximately the experimental p t distribution of pions and protons, the transverse 
momentum of the cluster (A(1900) in this case) has to be of the order of 2 GeV/c. Such large initial p t causes 
the final particles to be close in r] and a(r]) is only about 0.6. For other resonance cascades with larger values of 
AM in intermediate decays, the initial pt, cluster is smaller and thus a(r]) is larger. 

The two cluster models were used to generate several samples of 100,000 events containing particles from many 
clusters. The multiplicity of generated charged particles was counted in appropriate 77 bins and values of a 2 (C) 
were calculated. These are compared with the data from the analysis [7] of Au+Au collisions at ^/s NN = 200 GeV. 
In Fig. 3 (left plot) we show the dependence of cr 2 (C) on A??, the width of the -q bins. The ICM results are 
presented for two different values of effective cluster size, fc e ff=3.4 and £; e fr=4.3, which fit the experimental results 
for central (0-20%) and semi-peripheral (40-60%) Au+Au collisions, respectively. The same cluster sizes are then 
used to obtain ICM predictions for the dependence of a 2 (C) on the position in r\ of the bins, all with a width 
of 0.5 (Fig. 3, right plot). In the case of the RCM, only one value of fc e ff=3 is used and results are presented 
in Fig. 3. The dependence of cr 2 (C) on A-q shows that in peripheral collisions the effective cluster size is larger 
by about one. This fc e ff difference may be smaller if the clusters in central collisions are narrower in 77, as the 
dependence of er 2 (C) on 77 suggests. The two models (but especially ICM) described in this study represent the 
extreme cases of the narrowest and the widest possible pseudorapidity clusters. The results allow us to set lower 
(from RCM) and upper (from ICM) limits for the cluster size (according to Fig. 3, left plot, they are ~2.5 and 
^4.3 respectively). In a more realistic model, a mixture of clusters with different widths and sizes should be 
used, leading to k c ff values between the RCM and ICM predictions. However, the data from forward-backward 
correlations are not precise enough to determine fc e ff and cluster 77-width at the same time. 

5 Conclusions 

A large sample of the most central 3% of Au+Au collisions at y/s^ = 200 GeV was examined in the search for 
events with an unusual shape of the dN/drj distribution. After pileup removal only a few such events were left, 
which allows us to set an upper limit for the probability of such fluctuations at 10~ 5 . 

The forward-backward multiplicity correlations are sensitive to multiplicity fluctuations at short distances in 
77, which are expected if particles are produced in clusters. Using two simple clusters models, we found that the 
effective cluster size consistent with the experimental data is between 2.5 and 4. 
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